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Abstract 

Next-to-leading-order QCD analyses of the ZEUS data on deep inelastic scatter- 
ing together with fixed-target data have been performed, from which the gluon 
and quark densities of the proton and the value of the strong coupling constant, 
as{Mz), were extracted. The study includes a full treatment of the experimen- 
tal systematic uncertainties including point-to-point correlations. The resulting 
uncertainties in the parton density functions are presented. A combined fit for 
as{Mz) and the gluon and quark densities yields a value for as{Mz) in agree- 
ment with the world average. The parton density functions derived from ZEUS 
data alone indicate the importance of HERA data in determining the sea quark 
and gluon distributions at low x. The limits of applicability of the theoretical 
formalism have been explored by comparing the fit predictions to ZEUS data at 
very low Q^. 
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1 Introduction 



Studies of inclusive differential cross sections and structure functions, as measured in 
deep inelastic scattering (DIS) of leptons from hadron targets, played a crucial role in 
establishing the theory of perturbative quantum chromodynamics (pQCD). The next-to- 
leading-order (NLO) DGLAP evolution equations [1,2,3,4] form the basis for a successful 
description of the data over a broad kinematic range. Thus parton distribution functions 
(PDFs) and the value of the strong coupling constant, as{Mz), can be determined within 
this formalism. The availability of HERA data has greatly increased the kinematic range 
over which such studies can be made. 

The MRST [5] and CTEQ [6] groups have used the most recent HERA data [7, 8] in 
global fits to determine PDFs and as{Mz)- In recent years, estimating the uncertainties 
on PDFs from experimental sources, as well as from model assumptions, has become an 
issue [9, 10, 11, 12, 13, 14, 15, 16]. The CTEQ group has made a detailed study of the 
uncertainties on the PDFs due to experimental sources, whereas MRST provide four sets 
of PDFs from fits done with different theoretical assumptions. The best fits of these 
groups differ somewhat, reflecting differences of approach. The HI collaboration has also 
considered the uncertainties on the gluon distribution and as{Mz) resulting from a fit to 
HI and BCDMS data [8]. 

In this paper, the ZEUS data from 1996 and 1997 [7] have been used, together with fixed- 
target data, to extract gluon and sea densities with much improved precision compared to 
earlier work that used the ZEUS 1994/95 data [17,10]. The fixed-target data are important 
for a precise determination of the valence distributions. All parton distributions have been 
extracted taking into account the point-to-point correlated systematic uncertainties of the 
input data. 

The value of as{Mz) was set to the world-average value, as{Mz) = 0.118 [18], for the 
determination of parton distributions in the standard fit (called ZEUS-S). The increased 
precision of the data also allows a determination of the value of as{Mz)- The correlations 
between the shape of the parton distribution functions and the value of as{Mz) have been 
fully taken into account by making a simultaneous fit to determine the parton distribution 
parameters and as{Mz)- This fit is called ZEUS-a^. 

One of the main topics of this paper is an evaluation of the experimental uncertainties on 
the extracted parton distribution functions and on the value of as{Mz)- The treatment 
of point-to-point correlated systematic uncertainties reflects knowledge that such uncer- 
tainties are not always Gaussian distributed. Model and theoretical uncertainties have 
also been estimated. 

The role of ZEUS data has been explored by making a fit using ZEUS data alone. The 
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ZEUS charged current e+p data from 1994-97 [19], and the charged and neutral current 
e~p data from the 1998/99 runs [20,21] were used, together with the 1996/97 e~^p neutral 
current data, to make an extraction of the PDFs independently of other experiments. 
This fit is called ZEUS-0. 

The extent to which the NLO DGLAP formalism continues to provide a successful de- 
scription of the data over an increased kinematic range was investigated by comparing 
the ZEUS-S fit to the ZEUS high very low [22]. The combination 

of the improved fit analysis and the increased precision of these data, compared to those 
used [23] in the previous study [17], allows a low-Q^ limit to be put on the applicability 
of the NLO DGLAP description of DIS data. 

The paper is organised as follows. In Section ^, some theoretical background is given. In 
Section 0, the NLO QCD fits to the ZEUS data and fixed-target DIS data are described, 
paying particular attention to the treatment of experimental uncertainties. In Section 
the standard ZEUS-S fit is compared to data and the extracted parton distribution func- 
tions including their experimental uncertainties are presented. The analysis is extended to 
evaluate as{Mz) in the ZEUS-a^ fit and uncertainties from experimental and theoretical 
sources are discussed. In Section |, parton densities from the ZEUS-0 fit are presented 
and, in Section ^ the limitations of the NLO DGLAP formalism are considered. Sec- 
tion 1^ contains a summary and conclusions. In the Appendix, various ways of treating 
systematic uncertainties are discussed and compared. 



2 Theoretical Perspective 

The differential cross section for neutral current (NC) e~^p deep inelastic scattering (DIS) 
is given in terms of the structure functions F2, and xF^ by 

The kinematic variables are Bjorken's x = Q'^/{2p ■ q) and the negative invariant-mass 
squared of the exchanged virtual boson, = — g^, where p is the four- vector of the target 
proton and q = k — k' is the difference of the four-vectors of the incoming and outgoing 
leptons. The variable y is defined hj y = {p ■ q)/{p ■ k) and Kt = 1 ± (1 — y)^. It is also 
useful to define W'^, the virtual boson-proton squared centre-of-mass energy, given by 

W^ = {p + qf = Q^^^^ + ml, 

where is the proton mass. The reduced cross section is defined as 

(Pa{e~^p) Q^x 



a 
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so that it is equal to F2 when and xF^ are neghgible. For values much below 
the Z°-mass squared, the parity- violating structure function, xF^, is negligible, since the 
cross sections are dominated by virtual photon exchange. Then, provided that W"^ is large 
enough that target-mass and higher-twist contributions may be neglected, the structure 
function F2 can be simply interpreted in LO QCD as the sum of the quark distribution 
functions weighted with the quark-charges squared. To the same approximation, Fl, the 
longitudinal structure function, is zero. At NLO, these structure functions are related 
to the parton distributions in the proton through convolution with the QCD coefficient 
functions. Since the ZEUS data extend to high Q^, the coefficient functions also include 
exchange [24] . Measurement of the structure functions as a function of x and yields 
information on the shape of the parton distributions and, through their dependence, 
on the value of the strong coupling constant, as{Mz)- 

Before HERA data became available, leading-twist perturbative expansions of QCD, as 
formulated in the DGLAP evolution equations, were found to describe fixed-target data 
adequately down to ~ 4 GeV^ and x ~ 10"^. The QCD evolution was typically started 
from Qo ~ 4 GeV^ or higher. Convenient functional forms of the parton distribution 
functions (PDFs) were input, at Qo, and fitted to the data. At small x, these were 
xf{x) ^ Ax^. The fits gave 5 ~ 0.5 for valence distributions and 5 ~ (fiat) for the sea 
and the gluon distributions. 

It was shown [25] as early as 1974 that, according to QCD, this behaviour cannot persist 
to infinitely small values of x. At some point, a much steeper rise of the gluon distribution 
is expected, leading to a steeply rising behaviour of F2. However, it was unclear at what 
X values such behaviour should begin. Hence, prior to HERA operation, most predictions 
specified a gentle rise at small x, as expected from soft Regge-like behaviour. The dramatic 
rise in F2 observed in the early HERA data [26,27], at x ~ 10"'^, ~ 15 GeV^ was 
therefore a surprise. Furthermore, later HERA data [28,29] showed that this rise persisted 
down to surprisingly low Q^, of order ~ 1.5 GeV^, where x ~ 5 x 10~^. Applications 
of pQCD using the NLO DGLAP formalism to data in this kinematic region have been 
reasonably successful, although there are several issues that could limit the applicability 
of this formalism. 

One question is whether only a few terms in the perturbative expansion are adequate, 
given the large values of as at low . Recent NLO DGLAP fits including HERA data 
have used starting values as low as ~ 1 GeV^. These fits have sea input distributions 
with 5 ~ —0.2. However, such fits require the gluon input distributions to be valence- 
like [30,31], or even negative [5], at small x. This calls into question the applicability of 
the DGLAP formalism at these low values of Q^. 

Furthermore, at the low x values accessed at HERA, large ln(l/x) terms, which are not 
included in the DGLAP formalism, could be important. If so, the treatment may need to 



3 



be amended by consideration of BFKL dynamics [32,33,34,35,36,37]. 

Finally, the high gluon density observed at higher could lead to gluons screening each 
other from the virtual-boson probe, requiring non-linear terms in the evolution equations. 
These act oppositely to the linear terms, such that gluon evolution is slowed down and 
may even saturate [38]. 

It is unclear where any of these effects become important. Presently, the range of appli- 
cability of the NLO QCD expansion is a matter to be resolved by experiment. To draw 
firm conclusions requires precision data and a careful analysis of the uncertainties on the 
predicted shapes of the parton distributions. In the present paper, the high-precision data 
from the ZEUS experiment and all fixed-target experiments for which full information on 
correlated systematic uncertainties is available have been used to extract the PDFs and 
a^iMz) and to investigate the range of applicability of the NLO DGLAP formalism. 

3 Description of NLO QCD fits 

This section gives the specifications of the ZEUS-S and ZEUS-a^ global NLO QCD fits 
to the new ZEUS cross-section data [7] and fixed-target DIS data. 

The fixed-target data were included to constrain the fits at high x and provide information 
on the valence distributions and the flavour composition of the sea. All high-precision 
fixed-target data sets for which full information on the correlated systematic uncertainties 
is available have been used: 

• F2 data on /i - p scattering from BCDMS [39], NMC [40], and E665 [41]; 

• deuterium-target data from NMC [40] and E665 [41]. These were included in order to 
have M, d flavour separation; 

• NMC data on the ratio F2 / F2 [42]. These determine the ratio of the d io u valence 
shapes; 

• the CCFR [43] xF^ data, from (anti-)neutrino interactions on an iron target. These 
give the strongest constraint on high-x valence PDFs. They are used only in the x 
range 0.1 < x < 0.65 in order to minimize dependence on the heavy-target corrections. 
The latter were performed according to the prescription of MRST [31]. These XF3 data 
are unaffected by the recent re-analysis of CCFR F2 data [44, 45] . 

The deuterium data were corrected to represent {n + p)/2 by the prescription of Gomez 
et al. [46]. The fit results were found to be insensitive to the specific prescriptions used 
for heavy-target and deuterium corrections. 
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The fits were performed at leading twist. Tlie following cuts were made on the ZEUS and 
the fixed-target data: 

• > 2.5 GeV^ was required to remain in the kinematic region where perturbative 
QCD is expected to be applicable; 

• W'^ > 20 GeV^ was required to reduce the sensitivity to target-mass [47] and higher- 
twist [48] contributions which become important at low W^. 

The kinematic range covered by the data input to the fits is 6.3 x 10~^ < a; < 0.65 and 
2.5 < g2 < 30,000 GeV^. 

The QCD predictions for the structure functions needed to construct the reduced cross 
section were obtained by solving the DGLAP evolution equations at NLO in the MS 
scheme [49,50,51] with the renormalisation and factorisation scales chosen to be Q^. 
The DGLAP equations yield the quark and gluon momentum distributions (and thus 
the structure functions) at all values of Q^, provided they are input as functions of x at 
some input scale, Q^. The input scale was chosen to he Ql = 7 GeV^; however, there 
is no particular significance to this value since backward evolution was performed to fit 
lower-Q^ data. The choices of the value of Q^, the forms of the parameterisations of the 
parton distributions at Q^, and the cuts on the data to be fitted have all been varied in 
the course of systematic studies (see Section 4.4). 



3.1 Parameterisation of parton distribution functions 

The parton distribution functions (PDFs) were parameterised at by the form 

xf{x) = pix^^{l — xy^{l + p^x) 

so that the distributions are either zero, or singular, as x — 0, and tend to zero as 
X —* 1. The parton momentum distributions that were parameterised are: u valence, 
xu^{x)] d valence, xd^^x); total sea, xS{x); gluon, xg{x)] and the difference between the 
d and u contributions to the sea, xA = x{d — u). The total sea at Ql is made from the 
flavours up, xUsca.{x), down, xdsca.{x), strange, xSgeai^) and charm, xcgcai^), as follows 

XSsca.{x) = 0.2xS{x) 

OAxS^x) — 0.5xCsca(a;) — xA{x) 

OAxS{x) — 0.5xcsca{x) + xA{x) 

where the symbols Usea., c^sea, Ssca, Csea include both quark and antiquark contributions to 
the sea for each flavour. The charmed sea is generated as described in the next section. 



■^^sca ('^) 
•^dgea (•^) 



5 



The suppression of the strange sea to 20% of the total sea is consistent with neutrino- 
induced dimuon data from CCFR [52]. The fit results are insensitive to this assumption. 

The following parameters were fixed: 

• pi for xu^ and xd^, were fixed through the number sum-rules and pi for xg was fixed 
through the momentum sum-rule; 

• P2 = 0.5 was fixed for both valence distributions, since, after the cut x > 0.1 on 
the xF^ data little information on the low-x valence shapes survives. Allowing these 
parameters to vary, and varying the value of the low-x cut, produces values consistent 
with 0.5 and has negligible effect on the shapes of these distributions; 

• the only free parameter for the xA distribution is its normalisation, pi, because there is 
insufficient information on its shape without including E866 Drell-Yan data [53]. Thus, 
P2(A) = 0.5, p3(A) = p3(Sea) + 2 were fixed, following MRST [54,55,31], andp5(A) = 
0; the normalisation pi(A) was found to be compatible with the measured value of 
the Gottfried sum-rule [56,57]. The fit results are insensitive to these assumptions; 

• for the gluon distribution, p^ was set to zero, since this choice constrains the high- a; 
gluon to be positive. Allowing this parameter to vary in the fit produces values which 
are consistent with zero. 

There are thus 11 free parameters in the ZEUS-S fit, when the strong coupling constant 
is fixed to as{Mz) = 0.118 [18], and 12 free parameters in the ZEUS-tts fit. In the 
DGLAP evolution equations at NLO, as(Q^) is calculated to 2-loop accuracy. The evo- 
lution was performed with the program QCDNUM [58]. The evolution equations were 
written in terms of quark-flavour singlet and non-singlet distributions (made from the 
sea and valence quark distributions) and the gluon momentum distribution. These must 
be convoluted with coefficient functions in order to calculate structure functions. The 
coefficient functions are specific to the heavy-quark formalism used, as discussed below. 

3.2 Treatment of heavy quarks 

The treatment of the heavy-quark sea needs careful consideration. Many early global 
fits [59, 60, 61, 54, 30, 62, 63, 64, 65] used zero-mass variable-fiavour-number (ZMVFN) 
schemes, where, for example, the charmed quark (of mass rric) is only produced once 
> Ami; larger Q^, the charm distribution is generated by the splitting g ^ cc using 
the equations for massless partons. This is incorrect at threshold. Other authors [66,67,68] 
have used a fixed-fiavour-number (FFN) scheme, in which a cc pair is created by boson- 
gluon fusion for W'^ > (2mc + rrip)'^, (a W'^ which may correspond to <^ 4m^, if x 
is small) but charm is then treated as a heavy quark which is dynamically generated for 
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all Q^. There is then no concept of a charmed parton distribution and thus \a{Q'^ /ml) 
terms remain in the NLO BGF coefficient functions, since they cannot be summed and 
absorbed into the definition of the charm distribution. This is incorrect at high Q^. 
Recently several groups [69,70,71,72,73,74,75,76] have tried to construct general-mass 
variable-fiavour-number schemes which behave correctly from threshold to large . In 
this analysis, the scheme of Thorne and Roberts [77, 78, 79, 80] (TRVFN) has been used 
to interpolate between correct threshold and correct large-Q^ behaviour. The results are 
compared to those obtained using the FFN and ZMVFN schemes in Section O. 



3.3 Definition of and treatment of correlated systematic un- 
certainties 

The minimisation and the calculation of the covariance matrices were based on MI- 
NUIT [81]. The definition of the 'was 

2 >^ [Fiip, g) - F,(meas)]^ , 2 (-^^ 
^ = 2^ — 7^^— + (1) 

where 

A 

The symbol Fj(meas) represents a measured data point (structure function or reduced 
cross section) and the symbols (Xj^stat and cTj^unc represent its error from statistical and 
uncorrelated systematic uncertainties, respectively. The symbol F^^'^^'^^ i^p) represents 
the prediction from NLO QCD in terms of the theoretical parameters p (PDF parameters 
and as{Mz))- This prediction is modified to include the effect of the correlated systematic 
uncertainties as shown in Eq. @). The one-standard-deviation systematic uncertainty 
on data point i due to source A is referred to as A^^*^ and the parameters s\ represent 
independent Gaussian random variables with zero mean and unit variance for each source 
of systematic uncertainty. These parameters sx were fixed to zero to obtain the central 
values of the theoretical parameters, but they were allowed to vary for the error analysis, 
such that in addition to the usual Hessian matrix, Mj^, given by 



2 dpjdpk ' 

which is evaluated with respect to the theoretical parameters, a second Hessian matrix, 
CjA, given by 

^ 2 dpjdsx 
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was evaluated. The systematic covariance matrix is then given by = M~^CC'^M~^ [82] 
and the total covariance matrix by V^°^ = + V^^, where = M~^. Then the uncer- 
tainty on any distribution F may be calculated from 



by substituting V^, V^"^ or V^°^ for V , to obtain the statistical (and uncorrelated system- 
atic), correlated systematic or total experimental error band, respectively. This method of 
accounting for systematic uncertainties is equivalent, to first order, to the 'offset method', 
in which each sx is varied by its assumed uncertainty (±1), a new fit is performed for each 
of these variations, and the resulting deviations of the theoretical parameters from their 
central values are added in quadrature [16]. Either of these methods of treating system- 
atic uncertainties results in more conservative error estimates than alternative methods 
discussed in the Appendix. 

The normalisations of the data sets were taken as published, apart from BCDMS data, 
which were scaled down [83,54,55,31,30] by 2%. However, the normalisation uncertainties 
were included among the correlated systematic uncertainties. In total, 71 independent 
sources of systematic uncertainty were included (see Table p. 

4 Fit results, theoretical and model uncertainties and 
the extraction of as{'M.z) 

4.1 Fit quality and fit predictions 

The ZEUS-S fit, with a,(Mz) = 0.118, is shown in Figs. |ll-|. In Fig. § the fit prediction 
for F2 is shown compared to the ZEUS and proton fixed-target data as a function of x 
at low Q^. In Fig. |^, this comparison is made as a function of for x values in the 
range 6.3 x 10~^ < x < 0.65. For the fixed-target data, only the 7-exchange process 
contributes to F2, whereas, at high-Q^, there are also contributions from Z° exchange 
and 7/2'° interference. Thus, for comparability with the fixed-target results, the ZEUS 
data shown in these figures represent only that part of F2 due to 7 exchange, as denoted 
by the symbol Fj^"^. The fit gives an excellent description of the data. 

The goodness of fit cannot be judged from the calculated from statistical and uncor- 
related systematic errors alone. Re-evaluating the for the parameters resulting from 
the ZEUS-S fit by adding the statistical, uncorrelated and correlated systematic errors 
in quadrature gives a total P^r data point of 0.95 for 1263 data points and 11 free 
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parameters. The P^r data point for individual data sets calculated in the same way 
are listed in Table |I[ 

In Fig. 1^, the fit is compared to the ZEUS high-Q^ neutral current e~^p data. This figure 
also shows predictions for the neutral current e~p data [21], which were not included in 
the fit. The charged current e~^p [19] and e~p [20] data (which were also not included in 
the fit) are compared to the fit prediction in Fig. ^. These high-Q^ data are very well 
described by the fit. 

4.2 Parton distribution functions and Fl 

The PDF parameters extracted from the ZEUS-S fit at Qq = 7 GeV^ are given in Table ^ 
and the corresponding parton distributions at = 10 GeV^ are shown in Fig. ^(a). The 
precision of these distributions is considerably improved in comparison to a fit [10] using 
earlier ZEUS data. The total error band is dominated by systematic uncertainties. In 
Fig. ^(b), the ZEUS parton distributions are compared to the latest distributions from 
MRST [5] and CTEQ [6]. The differences between these sets of parton distributions are 
compatible with the size of the error bands on the ZEUS parton distributions. 

The PDFs extracted from the ZEUS-S fit are now considered in more detail. In these 
distributions, the contribution to the error bands coming from variation of as{Mz) will be 
indicated in addition to the contributions of correlated and uncorrelated experimental un- 
certainties. This additional uncertainty has been taken into account with full correlations 
by allowing as{Mz) to be a parameter of the ZEUS-a^ fit (see Section ^73|) . 

The valence distributions xu^ and xdy extracted from the fit are shown for several dif- 
ferent values in Figs. |^ and |^. The abscissa is linear and the ordinate logarithmic to 
illustrate the high-x behaviour of these valence distributions, where they are constrained 
by the fixed-target data. The distributions for = I GeV^ are obtained by backward ex- 
trapolation. The uncertainty is shown beneath each distribution in terms of the fractional 
differences from the central value. The w-valence distribution is much better determined 
than the d-valence distribution, since structure-function data from fixed-target experi- 
ments are dominated by the u quark. 

The extracted sea distribution and its uncertainty is shown for several values in Fig. ||. 
The uncertainty in these distributions is less than ~ 5% for > 2.5 GeV^ and 10~^ < 
X < 10~^, but considerable uncertainty remains for x > 0.1. The sea distribution rises at 
small X, even at = 1 GeV^. 

The corresponding gluon distribution and its uncertainty is shown for several values 
in Fig. p. The general shape of the error bands, with a narrowing at x ~ 0.1, is a 
consequence of the momentum sum rule. The gluon distribution is determined to within 
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~ 10% for Q'^ > 20 GeV^ and 10~^ < x < 10"^; and its uncertainty decreases as Q'^ 
increases. Considerable uncertainty remains for x > 0.1. The distribution rises steeply 
at low X for ^ 5 GeV^; however, at lower Q^, the low-x gluon shape is flatter. When 
the fit is extrapolated back to = 1 GeV^, the shape becomes valence-like and tends to 
become negative at the lowest x, although remaining consistent with zero. 

The shapes of the gluon and the sea distributions are compared in Fig. [TO. For > 
5 GeV^, the gluon density becomes much larger than the sea density, but for lower 
the sea density continues to rise at low x, whereas the gluon density is suppressed. The 
present analysis shows this contrasting behaviour of the low-x, low-Q^ gluon and sea 
distributions even more clearly than the previous study of earlier ZEUS data [17]. 

It is also interesting to compare the behaviour of the gluon and the sea NLO densities as 
a function of for fixed x values. This is shown in Fig. |ll]. The scaling violation of the 
gluon distribution at small x is striking, reflecting the singular behaviour of the Pgg and 
Pgg splitting functions as x ^ 0. 

The tendency of the gluon distribution to become negative at low x and low could be a 
signal that NLO QCD is inadequate in this kinematic region. However, the only physical 
requirement is that the structure functions calculated from the parton distributions are 
positive. Thus it is important to investigate the fit prediction for F^, the structure function 
most closely related to the gluon [84] . This is shown in Fig. 0. It exhibits similar features 
to the gluon. This will be discussed further in Section 0. 



4.3 The extraction of as{Mz) 

In the evolution of singlet quark distributions at intermediate x (0.01 < x < 0.3), the value 
of as{Mz) and the gluon shape are strongly correlated through the DGLAP equations, 
since an increase in as{Mz) can be compensated by a harder gluon distribution. This has 
restricted the precision of determinations of as{Mz) from NLO DGLAP fits to DIS data 
in the past. However, at small x [x < 0.01) this correlation is weakened, since the gluon 
then drives the behaviour of F2 as well as that of dF2/d\n{Q'^). Thus, precision low- a: 
data can be used to make a simultaneous fit for as{Mz) and the PDF parameters. In the 
ZEUS-Og fit, as{Mz) was left free, leading to 

as{Mz) = 0.1166 ± 0.0008(uncorr.) ± 0.0032(corr.) ± 0.0036(norm.), (3) 

where the three uncertainties arise from the following: statistical and other uncorrelated 
sources; correlated systematic sources from all the contributing experiments except that 
from their normalisations; the contribution from the latter normalisations. 

The difference between this value of as{Mz) and the value 0.118 used in the ZEUS-S 
fit does not produce any significant shifts in the PDF parameters as compared to those 



10 



determined in the ZEUS-S fit. However, the correlation between as{Mz) and the PDF 
parameters does increase their experimental uncertainties, particularly that of the gluon, 
as illustrated in Figs. P-|T2. 



4.4 Model uncertainties 

Sources of model uncertainty within the theoretical framework of leading-twist NLO QCD 
are now considered. The sensitivity of the results to the variation of input assumptions 
has been quantified in terms of the resulting variation in as{Mz), since it is the most 
sensitive parameter. 

Table ^summarizes the effect of varying the value of Q^, and the minimum Q^, x and 
(Qmin; ^min, W^mj^) of data entering the ZEUS-a^ fit, in terms of the shift in the central 
value of as{Mz). These variations produce only a small model uncertainty in as{Mz) and 
in the PDF parameters. 

It is also necessary to consider varying the form of the input PDF parameterisations. 
Variation in the gluon parameterisation produces the most significant effects since it is 
least well known. Allowing the high-x gluon to take a more complex form, with ^ 0, 
resulted in a shift of Aas{Mz) = +0.0002. Extending the form of the parameterisation 
from (1 +P5x) to (1 +p4y/x + p^x) for both the gluon and the other parton distributions 
resulted in a shift of Aas{Mz) = +0.0008. Allowing p2 to be free for the valence distri- 
butions had no further effect on the value of as{Mz)- Finally, information from Tevatron 
high-i?r jet production [85,86] was used to constrain the high-x gluon [5]. The corre- 
sponding shift in the central value of as{Mz) was Aas{Mz) = +0.0006 and the shape of 
the gluon was shifted to be harder at high x. However, these shifts are well within the 
error estimates for both as{Mz) and the gluon PDF parameters. 

A further significant choice is that of the heavy-quark production scheme. Repeating the 
fit using the FFN scheme or the ZMVFN scheme produced shifts of Aa^(M^) = ±0.0010. 
Variation of the heavy-quark mass within the FFN and TRVFN schemes produced smaller 
shifts. The choice of the heavy-quark scheme also affects the shape of the gluon, such 
that the FFN scheme gluon is steeper at small x than the ZMVFN scheme gluon, with 
the TRVFN gluon in between. The size of these shifts is well within the error estimates 
of the gluon PDF parameters. 

Thus, the total model uncertainty on as{Mz) is considerably smaller than the errors from 
correlated systematic and normalisation uncertainties and leads to 

as{Mz) = 0.1166 ± 0.0008(uncorr.) ± 0.0032(corr.) ± 0.0036(norm.) ± 0.0018(model). 
The PDF parameters are much less sensitive to the model assumptions than is as{Mz)- 
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It follows that the error bands illustrated on the parton densities in Figs. ^[TT] represent 
reasonable estimates of the total uncertainties. 

4.5 Uncertainties in the theoretical framework 

While the uncertainty within the theoretical framework of leading-twist NLOQCD is 
rather well defined, it is much more difficult to decide on the uncertainty caused by 
reasonable variations in the framework. In this Section, two variations on the framework 
are estimated; the treatment of higher-twist terms and the choice of the renormalisation 
and factorisation scales, which gives an estimate of the importance of the higher-order 
terms in the pQCD expansion. 

The analysis was performed at leading twist and accordingly a hard W"^ cut was made 
to remove the region where higher- twist effects are known to be important. In order 
to evaluate if there are residual effects of higher twist at such large W"^, this cut was 
lowered to > 4 GeV^ and the SLAG data [87] were included^. A fit in which the 
leading-twist predictions for F2 were modified by a factor (l + Zij/Q^) was then performed, 
where hi, i = 1,10, are parameters determined in ten bins of x [88]. This modification 
was not intended to provide a thorough study of the higher-twist effects themselves, 
but only as an estimate of the uncertainty introduced by neglecting them. Hence, a 
simple form of the higher-twist contribution was used, in which xF^ was not modified and 
the higher-twist terms for deuterium and proton targets were assumed to be the same. 
The contribution of higher twist was found to be negligible for x < 0.005, small and 
negative for 0.005 < x < 0.5 and large and positive for x > 0.5, where target-mass effects 
are important. Having determined the hi parameters in this fit, these parameters were 
fixed and a fit was performed with the usual hard W"^ cut (excluding SLAG data). This 
produced a shift of Aa,(Mz) = -0.0032. 

Variation of the renormalisation and factorisation scales used in the fit was also con- 
sidered. The choice of for these scales is conventional in the inclusive DIS process, 
and their variation is used as a crude way of estimating the importance of higher-order 
terms [89,90,91,92]. These scales were varied from 2Q^, independently and 

simultaneously. This produced shifts Aas{Mz) ~ ±0.004, mostly from the change in 
renormalisation scale. The result of making larger scale changes, such as — > 4Q^, is 
not presented because such large scale changes produce fits with much larger x^, which 
are unacceptable according to the 'hypothesis testing' criterion (see Appendix). It is un- 
clear that such arbitrary scale changes give any reasonable estimate of the importance 

^ Note that the for these data must be calculated by adding statistical and systematic errors in 
quadrature, since information on correlated point-to-point systematic uncertainties is not available. 
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of higher-order terms [5]. Several groups [93,94,95] have compared NLO and approxi- 
mate NNLO analyses. The change in as{Mz) obtained in these studies is in the range 
-0.0035 < Aa,(Mz) < -0.001. 

The uncertainties discussed in this Section are rather large. However, since these investi- 
gations are far from exhaustive and given the difficulties in defining a reasonable variation 
in the theoretical framework, they are not included in the uncertainties quoted on the final 
value of as{Mz). 



5 Parton densities from ZEUS data alone 

The fit using ZEUS data only (ZEUS-0) uses the charged current e~^p data [19] and 
the neutral and charged current e~p data [21,20] in addition to the e~^p neutral current 
data [7]. These high-Q^ data are very well described by the ZEUS-S fit, as illustrated 
in Figs. and |. However, in the ZEUS-0 fit these additional data sets were used 
instead of the fixed-target data to constrain the valence distributions. Note that the 
for these additional data sets must be calculated by adding statistical and systematic 
errors in quadrature. The correlated point-to-point systematic uncertainties are small 
compared to the statistical uncertainties for these data sets. Since the exclusion of the 
fixed-target data leaves no constraint on the flavour content of the sea, the value of pi 
for the A distribution was fixed to the value determined in the ZEUS-S fit. The value 
as{Mz) = 0.118 was fixed; all other parameters were varied as usual. 

The gluon and the sea distributions extracted from the ZEUS-0 fit are shown in Fig. |T^. 



Comparing this figure to Fig. [TO, it is clear that the gluon and sea densities are mainly 
determined by the ZEUS data for x < 10~^. The ZEUS-0 fit gives almost as good a 
determination of these distributions as the ZEUS-S fit over most of the x, plane used 
in the fit. 

The valence distributions extracted from the ZEUS-0 fit are shown in Fig. |1^. They are 
determined to a precision about a factor of two worse than in the ZEUS-S fit. The u- 
valence distribution is well determined; however, the o?- valence distribution is much more 
poorly determined. In the ZEUS-0 fit, the d-valence distribution is determined by the 
high-Q^ e+p charged current data. In contrast in the ZEUS-S fit the rf- valence distribution 
is determined by the deuterium fixed-target data. Recently it has been suggested that such 
measurements are subject to significant uncertainty from deuteron binding corrections [96, 
97,98,99,100]. The ZEUS-0 extraction does not suffer this uncertainty. It produces a 
larger (i- valence distribution at high-x than the ZEUS-S fit, as can be seen by comparison 
with Fig. 1^, but there is no disagreement within the limited statistical precision of the 
current high-Q^ data. 
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6 The transition to very low 



The ZEUS-S and ZEUS-a, fits and the NLO QCD fits of MRST [30, 31, 55, 5] and 
CTEQ [64,65,6] give good descriptions of data down to values of 1 — 2 GeV^. 
For such fits to be valid, it is necessary to assume that the formalism is valid even for low 
[101], where is large and perturbation theory may break down, as well as for very 
low X, where ln(l/x) resummation terms should become important [90,91,92,102,103]. 
High-density and non-perturbative effects [104] are also neglected. To investigate if there 
is a low-Q^ limit to the applicability of the NLO QCD DGLAP formalism, the ZEUS-S 
fit was extrapolated into the region covered by ZEUS shifted- vertex (SVX) data [17] 
and the precise ZEUS beam-pipe-tracker (BPT) data [22]. 

In Fig. |l^, the ZEUS 96/97 data and the SVX and BPT data are shown at very low 
compared to the predictions of the ZEUS-S fit. The increased precision of the new 
data, both at intermediate and at very low Q^, lead to a firmer conclusion than in the 
previous study [17]. The ZEUS-S fit is able to describe the data down to = 1.5 GeV^, 
but exceeds the data at = 0.8 GeV^, and clearly fails for Q'^ < 0.65 GeV^, even when 
the conservative error bands on the fit due to the correlated systematic uncertainties are 
included. Thus, the NLO DGLAP formalism describes the extreme steepness of the ZEUS 
data at intermediate (2.7 <Q^ < 200 GeV^) but is unable to accommodate the rapid 
transition to a flatter behaviour at < 1 GeV^. The ZEUS-S fit predictions for for 
very low-Q^ values are also shown in Fig. |1^. The significantly negative values of F^ for 
< 1 GeV^ are a further indication that the NLO DGLAP formalism is not applicable. 

7 Summary and conclusions 

The NLO DGLAP QCD formalism has been used to fit the 1996-1997 ZEUS data and 
fixed-target data in the kinematic region Q'^ > 2.5 GeV^, 6.3 x 10~^ < x < 0.65 and W'^ > 
20 GeV^. Full account has been taken of correlated experimental systematic uncertainties. 
A good description of the structure function and reduced cross section over the range 
from 2.5 to 30, 000 GeV^ has been obtained. 

The parton distribution functions for the u and d valence quarks, the gluon and the total 
sea have been determined and the results are compatible with those of MRST2001 and 
CTEQ6. The ZEUS data are crucial in determining the gluon and the sea distributions 
and a fit to ZEUS data alone shows that these data also constrain the valence-quark 
distributions. The new high-precision data allow a greatly improved determination of the 
gluon and sea distributions. 

At ~ 1 GeV^, the fit predicts that the sea distribution is still rising at small x, 
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whereas the gluon distribution is suppressed. The fit is unable to describe the precise 
ZEUS BPT data for < 1 GeV^ and also predicts unphysical negative values for Fl in 
this region. Hence the use of the NLO QCD DGLAP formahsm at < 1 GeV^ is 
questionable. 

The ZEUS data at low x have been used to extract the value of as{Mz) in a simultaneous 
fit to as{Mz) and the shapes of the input parton distributions, including correlations 
between them, giving 

as{Mz) = 0.1166 ± 0.0008(uncorr.) ± 0.0032(corr.) ± 0.0036(norm.) ± 0.0018(model). 

Uncertainties in the leading twist NLO QCD framework are also significant but cannot 
be easily quantified. 

The statistical accuracy of the ZEUS data is now sufficient to give a very good determi- 
nation of the sea and gluon PDFs. With the full HERA-II data sample, it will be possible 
to extend this analysis to give an accurate determination of all proton PDFs within a 
single experiment. 

Appendix 

A. Comparison of different ways of calculating 

The used in this analysis is defined in Eq. (||) and the modification of the theoretical 
predictions to account for correlated systematic uncertainties is given in Eq. (H). The 
has been evaluated with the systematic-offset parameters set to zero, sx = 0, with 
the consequence that the fitted theoretical predictions are as close as possible to the 
central values of the published data. The offset parameters were then allowed to vary 
in the evaluation of the error to account for correlation between systematic-uncertainty 
parameters and theoretical parameters, as described in Section p.3| . 

This method is referred to as the 'offset method', since it is approximately equivalent to 
offsetting each systematic parameter s\ by ±1, performing a new fit for each of these 
variations and adding in quadrature the resulting deviations of the theoretical parameters 
from their central values [16]. This procedure does not assume that the systematic errors 
are Gaussian distributed. This is a conservative method of error estimation as compared 
to the Hessian methods described below [11,16]. 

An alternative procedure would be to allow the systematic uncertainty parameters sx to 
vary in the main fit when determining the values of the theoretical parameters. This 
was the procedure adopted by a recent HI analysis [8], in which only HI and BCDMS 
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data were considered. This method is referred to as 'Hessian method 1'. The errors 
on the theoretical parameters are calculated from the inverse of a single Hessian matrix 
which expresses the variation of with respect to both theoretical and systematic offset 
parameters. Effectively, the theoretical prediction is not fitted to the central values of the 
published experimental data, but allows these data points to move collectively according 
to their correlated systematic uncertainties. The theoretical prediction determines the 
optimal settings for correlated systematic shifts of experimental data points such that 
the most consistent fit to all data sets is obtained. Thus the fit correlates the systematic 
shifts in one experiment to those in another experiment. 

Hessian method 1 becomes a cumbersome procedure when the number of sources of 
systematic uncertainty is large, as in the present global DIS analysis. Recently the 
CTEQ [13, 14] collaboration has given an elegant analytical method for performing the 
minimization with respect to systematic-uncertainty parameters. This gives a new for- 
mulation of the x^: 



such that the uncorrelated and correlated systematic contributions to the can be eval- 
uated separately. This method is referred to as 'Hessian method 2'. 

These two Hessian methods have been compared for the ZEUS-0 fit, in which the sys- 
tematic uncertainties are well understood. The results are very similar, as expected if 
the systematic uncertainties are Gaussian and the values A^|* are standard deviations. 
However, if data sets from different experiments are used in the fit, the results of these 
two Hessian methods are only similar if normalisation uncertainties are not included. 

The offset method has been compared to Hessian method 2 by performing the ZEUS-a^ fit 
to global DIS data using Hessian method 2 to calculate the x^. Normalisation uncertain- 
ties were excluded and as{Mz) was included as one of the theoretical parameters. This 
fit yields as{Mz) = 0.1120 ± 0.0013, where the error represents the total experimental 
uncertainty from correlated and uncorrelated sources, excluding normalisation uncertain- 
ties. Thus this value should be compared with as{Mz) = 0.1166 ±0.0033, evaluated using 
the offset method, also excluding normalisation uncertainties (see Eq. Hessian method 




where 




and 
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2 gives a much reduced error estimate for both as{Mz) and the PDF parameters. The 
value of as{Mz) is shifted from that obtained by the offset method. The PDF parameters 
are not affected as strongly; their values are shifted by amounts which are well within the 
error estimates quoted for the offset method. 

To compare the of fi^s performed using the offset method and Hessian method 
2, it is necessary to use a common method of calculation. Table ^ presents the 
for the theoretical parameters obtained using each of these methods, re-evaluated by 
adding statistical and systematic errors in quadrature. For both methods, as{Mz) has 
been included among the theoretical parameters and normalisation uncertainties have 
not been included among the systematic parameters. The total increase of for Hessian 
method 2 as compared to the offset method is Ax^ = 283. Thus the results of Hessian 
method 2 represent a fit with an unacceptably large value of when judged in this 
conventional way. 



B. Parameter estimation and hypothesis testing 

To appreciate the signficance of the difference in between various fits, the distinction 
between the changes appropriate for parameter estimation and for hypothesis testing 
should be considered. Assuming that the experimental uncertainties that contribute have 
Gaussian distributions, errors on theoretical parameters that are fitted within a fixed the- 
oretical framework are derived from the criterion for 'parameter estimation' x^ —>■ Xmin + 1- 
However, the goodness of fit of a theoretical hypothesis is judged on the 'hypothesis test- 
ing' criterion, such that its x^ should be approximately in the range N± v^2A^), where 
is the number of degrees of freedom. 

Fitting DIS data for PDF parameters and as{Mz) is not a clean situation of either param- 
eter estimation or hypothesis testing, nor are the contributing experimental uncertainties 
always Gaussian distributed. Within the theoretical framework of leading-twist NLO 
QCD, many model inputs, such as the form of the PDF parameterisations, the values of 
cuts, the value of Ql, the data sets used in the fit, etc. can be varied. These represent 
different hypotheses and they are accepted provided the fit x^ falls within the hypothesis- 
testing criterion. The theoretical parameters obtained for these different model hypotheses 
can differ from those obtained in the standard fit by more than their errors as evaluated 
using the parameter-estimation criterion. In this case, the model error on the parameters 
can exceed the estimate of the total experimental error. This does not happen for the 
offset method, in which the uncorrelated experimental errors evaluated by the parameter- 
estimation criterion are augmented by the contribution of the correlated experimental 



systematic uncertainties, as explained in Section p.3i The shifts in theoretical parameter 



values for the different model hypotheses were found to be well within the total exper- 
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imental error estimates^. However, this is no longer the case when the fit is performed 
using Hessian method 2. 

The CTEQ collaboration [14, 6] have considered this problem. They consider that ~^ 
+ 1 is not a reasonable tolerance on a global fit to approximately 1200 data points 
from diverse sources, with theoretical and model uncertainties that are hard to quantify 
and experimental uncertainties that may not be Gaussian distributed. They have tried to 
formulate criteria for a more reasonable setting of the tolerance T, such that + 
becomes the variation on the basis of which errors on parameters are calculated. In 
setting this tolerance they have considered that all of the current world data sets must be 
acceptable and compatible at some level. The level of tolerance they suggest is T ~ 10. 
The error estimates of the present fit have been re-evaluated using Hessian method 2 
for various values of the tolerance. For T = 7 the errors on the PDF parameters and 
on as{Mz) are very similar to those of the offset method performed under the same 
conditions 0. For example, the result as{Mz) = 0.1120 ± 0.0033 is obtained. Note that 

the value T = 7 is similar to the hypothesis-testing tolerance T = \J~s/l2N) for the fits. 

Thus the offset method and the Hessian method with a modified tolerance T = 7 give 
similar error estimates. In choosing between these methods, there are some additional 
considerations. In the Hessian method 2, it is necessary to check that data points are 
not shifted far outside their uncertainties. When the ZEUS-a<j and the ZEUS-S fit are 
done by Hessian method 2 some of the systematic shifts for the ten classes of systematic 
uncertainty of the ZEUS data move by ~ ±1.4 standard deviations. No single kinematic 
region responsible for these shifts could be identified. Whereas these shifts are not very 
large, it is significant that they differ from the systematic shifts to ZEUS data determined 
in the CTEQ fit [6]. They also differ from those determined in the ZEUS-0 fit done by 
Hessian method 2. Making different model assumptions in the fits also produces somewhat 
different systematic shifts. It seems unreasonable to let variations in the model, or the 
choice of data included in the fit, change the best estimate of the central value of the data 
points. 

In summary, the offset method has been selected for several reasons. First, its fit results 
make theoretical predictions that are as close to the central values of the published data 
points as possible. The selection of data sets included in the fit or superficial changes 
to the model do not change the best estimate of the central value of the data points. 
Secondly, it is approximately equivalent to a method that does not assume that experi- 
mental systematic uncertainties are Gaussian distributed. Thirdly, its results produce an 
acceptable when re-evaluated conventionally by adding systematic and statistical er- 

^ Note that this is true whether or not normahsation uncertainties are included in these estimates. 
^ This remains the case when normalisation uncertainties are introduced into each of these methods. 
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rors in quadrature. Fourthly, its error estimates take account of the fact that the purpose 
is to estimate errors on the PDF parameters and as{Mz) within a general theoretical 
framework not specific to particular model choices. Quantitatively, the error estimates of 
the offset method correspond to those that would be obtained using the more generous 
tolerance T = 7 in the more statistically powerful Hessian methods. 
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Experiment 


Number of data 


per 


Number of correlated 


Number of normalisation 




points 


data point 


systematic uncertainties 


uncertainties 


ZEUS96/97 [7] 


242 


0.85 


10 


2 


BCDMS p [39] 


305 


0.94 


5 


5 


NMC p [40] 


218 


1.21 


12 


4 


NMC D [40] 


218 


0.92 




NMC D/p [42] 


129 


0.94 


5 





E665 D [41] 


47 


0.94 


7 


2 


E665 p [41] 


47 


1.16 




1 


CCFRxF3 [43] 


57 


0.40 


18 






Table 1: Table of for the data sets used in the ZEUS-S NLO QCD fit, evalu- 
ated by adding all systematic and statistical uncertainties in quadrature. Note that 
for CCFR data no separate total systematic uncertainty is supplied, so that this 
procedure overestimates the total uncertainty. The number of correlated systematic 
uncertainties for each data set is also given. Note that the systematic uncertainties 
for the p and D data sets of NMC and E665 must be taken together. The normal- 
isations of the four beam energies comprising the NMC data are the same for the 
p and D targets, whereas for E665 data there is a separate normalisation uncer- 
tainty for the p and D targets as well as a common normalisation uncertainty. The 
number of normalisation uncertainties for BCDMS data derives from the four beam 
energies of the data and an overall normalisation uncertainty. There are two ZEUS 
normalisation uncertainties: an overall uncertainty and the relative uncertainty of 
the data for which < 30 GeV^, with respect to the higher data. The CCFR 
normalisation uncertainty is included amongst its systematic uncertainties. 
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PDF 


Pi 


P2 


P3 


P5 




(1.69 ±0.01 ±0.06) 


0.5 


4.00 ±0.01 ±0.08 


5.04 ±0.09 ±0.64 




(0.96 ±0.01 ±0.08 


0.5 


5.33 ±0.09 ±0.48 


6.2 ±0.4 ±2.3 


xS 


0.603 ±0.007 ±0.048 


-0.235 ± 0.002 ±0.012 


8.9 ± 0.2 ± 1.2 


6.8 ±0.4 ±2.0 


xg 


(1.77 ±0.09 ±0.49) 


-0.20 ±0.01 ±0.04 


6.2 ± 0.2 ± 1.2 





xA 


0.27 ±0.01 ±0.06 


0.5 


(10.9 ± 0.2 ± 1.2) 






Table 2: Table of PDF parameters at Qq, as determined from the ZEUS-S fit. 
The first uncertainty given derives from statistical and other uncorrelated sources 
and the second uncertainty is the additional contribution from correlated systematic 
uncertainties. The numbers in parentheses were derived from the fitted parameters 
as described in the text. 
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Nominal value 


New value 




Ql = 7 GeV^ 


g2 = 4 GeV^ 
= 10 GeV^ 


+0.0008 
-0.0004 


= 2.5 GeV^ 


g^i„ = 4.5 GeV^ 


-0.0007 


a^min = 6.3 X 10"^ 


r ■ — 10~^ 


-0.0005 


W"^- =20 GeV^ 


W^- =10 GeV^ 

' mm " ^ 


-0.0005 



Table 3: Shifts in the central value ofas{Mz) with variation of the nominal input 
values for the ZEUS-ag fit. 



Experiment 


Data points 


X^/data point 


X^/data point 






Hessian method 2 


Offset method 


ZEUS96/97 


242 


1.37 


0.83 


BCDMS p 


305 


0.95 


0.89 


NMC p 


218 


1.50 


1.26 


NMC D 


218 


1.15 


0.96 


NMC D/p 


129 


0.97 


0.93 


E665 D 


47 


0.97 


0.94 


E665 p 


47 


1.17 


1.16 


CCFRxFS 


57 


0.99 


0.39 



Table 4: Table of calculated by adding systematic and statistical errors in 
quadrature for the theoretical parameters determined by the offset method and Hes- 
sian method 2 
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ZEUS 




Figure 1: The ZEUS-S NLO QCD fit compared to ZEUS 96/97 and proton 
fixed-target F2 data. The error bands of the fit represent the total experimental 
uncertainty from both correlated and uncorrelated sources. 
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Figure 2: The ZEUS-S NLO QCD fit compared to ZEUS 96/97 and proton 
fixed-target F2 data. The error bands are defined in the caption to Fig. |7|. 
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ZEUS 




Figure 3: The ZEUS-S NLO QCD fit compared to the ZEUS high-Q^ e+p and e'p 
neutral current reduced cross sections. The error bands are defined in the caption 
to Fig. Note that the e^p data were taken at y/s = 300 GeV, whereas the e"p 
data were taken = 318 GeV. 
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ZEUS 



0.75 
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0.25 



Figure 4: The ZEUS-S NLO QCD fit compared to the ZEUS high-Q'^ e+p and e'p 
charged current reduced cross sections. The error bands are defined in the caption 
to Fig. [|. Note that the e~^p data were taken at ^ s = 300 GeV, whereas the e~p 
data were taken at y/s = 318 GeV. 
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Figure 5: (a) The gluon, sea, u and d valence distributions extracted from the 
standard ZEUS-S NLO QCD fit at Q"^ = 10 GeV^ The error bands in this figure 
show the uncertainty from statistical and other uncorrected sources separately from 
the total uncertainty including correlated systematic uncertainties, (b) The gluon, 
sea, u and d valence distributions extracted from the ZEUS-S NLO QCD fit at = 
10 GeV^, compared to those extracted from the fits MRST2001 [5] and CTEQ6 [6]. 
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ZEUS 




Figure 6: The xu^ distribution from the ZEUS-S NLO QCD fit. The cross- 
hatched error bands show the statistical and uncorrelated systematic uncertainty, 
the grey error bands show the total experimental uncertainty including correlated 
systematic uncertainties (both evaluated from the ZEUS-S fit) and the hatched error 
bands show the additional uncertainty coming from variation of the strong coupling 
constant as{Mz) (evaluated from the ZEUS-ag fit). The uncertainties on these 
distributions are shown beneath each distribution as fractional differences from the 
central value. 
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ZEUS 




Figure 7: The xd^ distribution from the ZEUS-S NLO QCD fit. The error bands 
are defined in the caption to Fig. The uncertainties on these distributions are 
shown beneath each distribution as fractional differences from the central value. 
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ZEUS 




Figure 8: The sea distribution from the ZEUS-S NLO QCD fit for various 
values. The error bands are defined in the caption to Fig. The uncertainties 
on these distributions are shown beneath each distribution as fractional differences 
from the central value. 
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ZEUS 




Figure 9: The gluon distribution from the ZEUS-S NLO QCD fit for various 
values. The error bands are defined in the caption to Fig. |^. The uncertainties 
on these distributions are shown beneath each distribution as fractional differences 
from the central value. Note that this uncertainty is not shown when the central 
value of the gluon distribution becomes negative. 
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ZEUS 




Figure 10: Comparison of the gluon and sea distributions from the ZEUS-S NLO 
QCD fit for various values. The error bands are defined in the caption to Fig. 
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ZEUS 




Figure 11: (a) The gluon distribution from the ZEUS-S NLO QCD fit as a 
function of for fixed x values (h) The sea distribution from the ZEUS-S NLO 
QCD fit as a function of for fixed x values. The error bands are defined in the 
caption to Fig. ||. 
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ure 12: The longitudinal structure function Fl from the ZEUS-S NLO QCD 
The error hands are defined in the caption to Fig. 
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Figure 13: The gluon and sea distributions from the ZEUS-0 NLO QCD fit 
in various bins. The error bands show the uncertainty from statistical and 
other uncorrected sources separately from the total uncertainty including correlated 
systematic uncertainties. The value of as{Mz) = 0.118 is fixed. 
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ZEUS 




Figure 14: The xu^, and xd^ distributions from the ZEUS-0 NLO QCD fit. The 
error hands are defined in the caption to Fig. |7^. The value of as{Mz) = 0.118 is 
fixed. 
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Figure 15: F2 data at very low (including SVX95 and BPT97 data) compared 
to the backward extrapolated ZEUS-S NLO QCD fit. The error bands are defined 
in the caption to Fig. |i|. 



42 



ZEUS 




Figure 16: The predictions for at very low from the backward extrapolated 
ZEUS-S NLO QCD fit. The error bands are defined in the caption to Fig. 
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